The pairing Hamiltonian for one pair of identical nucleons bound in a potential well 
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The problem of one pair of identical nucleons sitting in j\f single particle levels of a potential 
well and interacting through the pairing force is treated introducing even Grassmann variables. 
The eigenvectors are analytically expressed solely in terms of these with coefficients fixed by the 
eigenvalues and the single particle energies. When the latter are those of an harmonic oscillator 
well an accurate expression is derived for both the collective eigenvalue and for those trapped in 
between the single particle levels, for any strength of the pairing interaction and any number of levels. 
Notably the trapped solutions are labelled through an index upon which they depend parabolically. 
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We have recently obtained, in the framework of the Grassmann algebra, the analytic expressions of the eigenvalues 
and the eigenvectors of a system of n pairs of like-nucleons interacting through the pairing Hamiltonian and sitting 
in one single-particle level |lj. 

We extend the analysis to the case of Af single-particle levels, with energies e\, e%, ■ ■ ■ ejv and angular momenta 
ji , j2, • • • jj\f (all the j's being assumed to be different): here the Hamiltonian, for identical particles, reads 
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the even (commuting, nilpotent) Grassmann variables. The latter is associated with a pair of fermions with vanishing 
third component of the total angular momentum (M—0). To start with we confine ourselves to consider one pair only. 

Notwithstanding the presence of both the A's and the 92's, the Hamiltonian ([l]) is diagonalized in the 2Af dimensional 
basis 
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which extends the one we have introduced in ref. and describes the two nucleons in the same single-particle level, 
fn (|^) 2£!„ = 2j v + 1 is the degeneracy of the level j v . 

For one pair, only states with seniority v=0 and 2 are allowed. In the basis (||) the eigenvalues of the v=2 states 
are trivial, whereas for those of the v=0 states one recovers the well-known secular equation 
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The corresponding eigenvectors are 



and 
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If a single particle level with j„ = 1/2 (fi„ = 1) is present, the dimension of the basis is actually 2j\f — 1, because, obviously, 
is identically zero. 
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with v — 1, • • - A/", the coefficients 
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fulfilling the system of equations 
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The above system is easily solved and yields the noticeable formula (see also 
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which shows that for a given set of single-particle energies the v=0 eigenvectors are fixed by the corresponding 
eigenvalues. When G is large (||) develops a collective solution with large E : hence the associated tend to become 
all equal and correspond to a coherent superposition of the so-called s-quasibosons. On the other hand, in the limit 
G — > 0, where E ~ 2e„, only one component of the basis, i.e. the z/-th one, survives in the wavefunction of the 
"trapped" states. 

In general the eigenvalues (and hence the eigenvectors) of (Q) stem from an interplay between the single-particle 
energies and degeneracies. Of course this interplay can be numerically explored. Here we pursue the scope analytically 
when the e v and the Q. v are available. One of the few cases where this occurs is for the harmonic oscillator well, where 
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Accordingly the secular equation (Q), for M levels, becomes 
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where G = G/2fiu and the energies are measured in units of %u> (2e}v = 2 AT + 3). 
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FIG. 1. The figure shows the solutions E of Eq.(|l(]|), for Af—6 (upper curves) and 8 (lower curves), as functions of G. One 
can see that with the harmonic oscillator well each trapped solution En for G > 0.1 tends approximatively to the single particle 
energy 2ejv. 

In Fig.l the numerical solutions of (|l^) are displayed for A/=6 and 8 versus G. Remarkably, the dependence upon 
G is seen to be lost for G > 0.1. Furthermore in this regime the trapped solutions are mildly depending upon by Af. 

Although aware that the solutions of ([Io|) cannot, in general, be algebraically expressed (for Af > 5), we explore 
whether the simple ansatz 
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provides a good representation of the trapped solutions (the collective solution E c will be separately treated). 
To fix the coefficients a, b and c, we recast (10) in the polynomial form 
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finding for the first three coefficients the expressions 
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Then the first three Viete equations, namely 
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yield a non-linear system in the unknowns a, b and c, if E c is known. This system can be solved by expressing, via 
eq.(|l^), c as a function of a and b 
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In turn (|l9j), inserted into (|l7|), yields b as a function of a. One finds 
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Finally, from (|18|), a cumbersome equation for a follows, not reported here. While (|16|) and ([L7|) are easily solved 
analytically, the non linear equation ( |18| ) for a can only be solved numerically and admits, in general, several solutions. 
The one appropriate for our problem is selected by requiring that the trapped energies lye in between the single-particle 
levels of the harmonic oscillators. Moreover, it should be pointed out that, owing to the high degree of non-linearity 
of (|l^), this solution turns out to be extremely sensitive to the collective energy E c , especially when G is large. 

Hence an accurate expression for the collective energy is needed. We look for the latter in the domain of small G 
(say G=0.05 - 0.1), as it follows from the empirical determination of G in atomic nuclei and from the experimental 
nuclear single-particle levels M. 
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For this scope we start by realizing that E c (0,Af)—3 and E c (Go,Af)=0, being 
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The above constraints are fulfilled by the cubic function (in G) 
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which thus provides an excellent representation of the collective energy (see Table 1). If an even better E c is wished, 
one can proceed perturbatively setting 



E c = E^ + S 

and expanding in the very small parameter S/M(N) where 

M(N) = 2N + 3- E { ° ] 

One thus gets the recursive relation 
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The energies provided by ( |28| ) fastly converge to the exact solution, as shown in Table 1, but, as above mentioned, a 
high precision is required, which obtains after 4 iterations. 
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-1.66966 


-1.47625 


-1.48025 
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-1.48025 


-1.48025 
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-8.37559 


-4.88888 


-5.44566 




-5.44568 


-5.44568 


-5.44568 
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-24.4931 


-3.26859 


-10.7011 




-11.0491 


-11.0546 


-11.0546 



TABLE I. Comparison between the exact E^ and the approximate E^ [eq. and eq.([28|)] collective energies for some 
values of Af and G=0.05 and 0.1 
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Formula ( |2§| ) yields E c also when G is large. Here however the analogous of (|25j) follows by expanding Eq.([L0|) in 
the parameter (2N + 3)/E. One thus gets for collective energy, expanded up to terms 1/G, the expression 

^.-f^ +1)( ^ +2) + ^ +1) - 2 ^)(^)^ . 

The above, when inserted in (^8|), yields results as accurate as those obtained in the domain of small G. 

With the collective energy fixed, the coefficients a, b and c can be found. We quote in Table 2, as an example, our 
predictions for the eigenvalues of the pairing hamiltonian for one pair in the J\f = 5 case, using as input ( p9| ) when 
G = 1 and G = 5 and J25| ) when G = 0.05 and G = 0.1. Our results are seen to agree with the exact ones obtained 
via the numerical solution of ( |l0| ) to better than 0.27%. This occurs as well for all the cases we have explored. Thus 
the simple ansatz ([ll]) is remarkably accurate. Furthermore the solutions (|ll"|), when G is large, scale in G. Indeed, 
in this condition, the right-hand sides of the Vieta equations ( p^[ - |l8|) are easily seen to be G-independent when the 
collective solution is given by (p9[). 
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TABLE II. Comparison between exact (e) "trapped" 


solutions of Eq.(|l0[) and approximate (app) 


ones, obtained from the 



ansatz ([Ll]) for J\f = 5 levels. The coefficients (a, b, c) of the parabola are (0.086,1.738,4.281) when G = 0.05, (0.014,2.183,3.427) 
when G = 0.1, (0.027,2.175,3.160) when G = 1 and (0.029,2.167,3.151) when G = 5. 



In conclusion, although the pairing problem can be, of course, solved numerically, yet we believe that our semi- 
analytical solution might be of some help for treating the situation when n pairs, sitting in an harmonic oscillator 
well, are present. Also interesting appears the extension of the present analysis to the situation where the pair is 
made out of a neutron and a proton, particularly when these are in an isospin singlet state |6) . In this case indeed the 
partners, in order to feel the pairing interaction, are forced to seat in different shells, at least 2huj apart. Whether in 
these conditions our semianalytical solution holds valid as well and a collective mode eventually develops is an issue 
we are currently exploring. 
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